Quantum state transfer via temporal kicking of information 
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We propose a strategy for perfect state transfer in spin chains based on the use of an unmodulated 
coupling Hamiltonian whose coefficients are explicitly time dependent. We show that, if specific and 
non-demanding conditions are satisfied by the temporal behavior of the coupling strengths, our 
model allows perfect state transfer. The paradigma put forward by our proposal holds the promises 
to set an alternative standard to the use of clever encoding and coupling-strength engineering for 
perfect state transfer. 

PACS numbers: 03.67.-a, 03.67.Hk, 75.10.Pq 



Multiple-spin systems, in particular spin chains, have 
recently been the object of extensive studies. From the 
quantum information processing (QIP) viewpoint, such 
systems embody valuable media for quantum protocols. 
In fact, it has been found that specific forms of built- 
in and permanent intra-register couplings, such as those 
typical of spin-chain models, could be used for the pur- 
poses of quantum computation |l[ and communication 
in quantum networks. This second possibility is par- 
ticularly interesting. In fact, while photons are ideal 
candidates for the long-haul transmission of informa- 
tion among different local nodes 0, their use in hybrid 
architectures for quantum networks (where the nodes 
are embodied by matter-like systems) requires the use 
of quantum interfaces. Despite the impressive experi- 
mental successes reported in this area, the realization of 
an interface is usually accompanied by technical prob- 
lems and errors in transmission. On the other hand, for 
short-distance quantum communication, an alternative 
arrangement where information carriers and processors 
are embodied by physical systems having the same na- 
ture could be more advantageous. In the seminal work in 
Ref. Q, the idea of using spin chains as quantum wires 
has been proposed as a way to avoid the interfacing prob- 
lems mentioned above. The original idea has then been 
extended along various directions (for more details, see 
Ref. 4]). In particular. Refs. [5] showed that, by properly 
engineering the strength of the couplings in the chain, 
perfect state transfer could be achieved. A proposal to 
bypass the initialization of the medium by means of local 
operations and measurements on the extremal spins of 
the chain has been recently put forward @. A similar 
result can be obtained by properly encoding the state to 
transmit in two spins at one end of the chain Q. 

Here, an alternative strategy to the one proposed in 
Ref. [5| is presented. Instead of the pre-arrangement of 
the interaction strengths across the chain, we consider 
a uniform distribution of time-dependent couplings. We 
show that perfect transfer of information is achieved in 
our time-dependent architecture as well and that a for- 
mal mapping of the dynamics achieved in the two models 
is in order, thus letting our proposal emerge as an al- 
ternative paradigma for perfect quantum state transfer. 



We also study a scenario where the coupling strengths 
of the Hamiltonian are constant, while an external time- 
dependent magnetic field is used in order to provide the 
necessary temporal modulation required for information 
transmission. We believe that the introduction of a time- 
dependent term in the intra-chain couplings allows to 
bypass the intrinsic rigidity of protocols based on un- 
modulated and pre-engineered spin media. In fact, it is 
usually the case that spin-chain models are arranged so as 
to implement a specific communication or computational 
task and cannot be "recycled" and used for a different 
one. Moreover, should an experimentally-prepared pat- 
tern of coupling strengths be found not accurate enough, 
the whole medium should be discarded (unless a non- 
ideal performance of the quantum protocol could be tol- 
erated). Differently, by allowing the presence of time- 
dependent Hamiltonian terms, a more dynamical "ad- 
justment" of the medium performances would be possible 
by implementing a simple feedback loop: state transmis- 
sion can be tested using a given functional form of the 
time-dependent terms of the Hamiltonian and, according 
to the results, this can be tuned so as to converge to- 
ward better performances. A time-dependent scheme for 
quantum state transfer has been discussed by Lyakhov 
and Bruder in Ref. [8| . Building up on general result de- 
scribed in Ref. , they only let the first and last spin-pair 
of a chain to experience time-dependent couplings. Our 
scheme, on the other hand, is different as we either con- 
sider a full set of time-dependent interaction strengths or 
a fully homogeneous pattern of constant couplings with 
a time-dependent external magnetic field. 

The remainder of this article is organized as follows. 
In Sec. U we describe the model used in our proposal. 
In Sec. |IT] we shortly review the main features of the 
information-flux approach, which is the main tool used 
in our investigation, and adapt it to our time-dependent 
analysis. Sec. Mil illustrates the way we obtain perfect 
state transfer when the coupling strengths depend on 
the interaction time via sharply rising/ lowering pulses, 
while Sec. HVl is dedicated to the investigation of a much 
more realistic setting that considers finite rising times 
and pulse duration. In Sec. [V] the proposed model is 
exploited in order to create a genuine multi-partite en- 
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tangled state. Finally, Sec. I VII summarizes our results. 



I. THE MODEL 

The system we analyze is an open spin-chain of N el- 
ements, whose Hamiltonian reads 

N-l N 

n = Yl [Jx{t)XiX i+1 + J y (t)YiY i+1 ] + B (t)Zi- (i) 
1=1 t=l 

Here, J x (t) and J y (t) are the coupling strengths of the 
pairwise interaction between adjacent spins and B(t) is 
a magnetic field. In our notation, X, Y and Z denote 
respectively the x, y and z-Pauli matrix. For the sake of 
simplicity, we consider N as an odd number. However, 
all the results presented here can be straightforwardly 
adapted to the case of an even number of spins in the 
chain. Physical units are chosen throughout the article 
so that h = 1. It is important to note that the inter- 
spin couplings and the amplitude of the magnetic field 
are site independent. This is a feature that diffcrcnciates 
the system at hand from a few previous proposals for 
perfect state transfer available in the literature 0-0, ■ 
Although Eq. ([1]) is spatially unmodulated (in analogy 
with Q ) , the price to pay in order to achieve unit transfer 
fidelity for any length of the chain is the time dependence 
of the interaction strengths. We will show later in this 
article that such a request can also be relaxed: perfect 
state transfer can be achieved even if only the magnetic 
field is time-dependent while any other term is constant. 

Eq. (JXJ) does not preserve the number of spin- 
excitations as it does not commute with the total z- 
componcnt of the spin of the system. Frequently, when 
one faces the case of dynamics ruled by spin-preserving 
models, on the assumption of proper boundary condi- 
tions, it is convenient to diagonalize the coupling Hamil- 
tonian by means of a sequence of operation comprising 
Wigner- Jordan, Fourier, and Bogoliubov transforms [TTj | . 
In our case, due to the spin-non preserving nature of 
%, the mutual coupling among sectors of the Hilbert 
space labeled by different quantum spin numbers has 
to be considered. Rather than applying techniques for 
the exact diagonalization of Eq. ([T]), here we tackle the 
evolution of the system by means of an information-flux 
approach, which is specifically designed for multi-spin in- 
teractions [HI, [III]. Our method does not rely on the 
explicit analysis of the energy spectrum of the chain and 
allows us to gather an intuitive picture of the dynam- 
ics at hand. Such an approach has already proved its 
flexibility and efficiency in identifying Hamiltonian con- 
figurations suitable for state transfer and entanglement 
generation @, 



II. INFORMATION-FLUX APPROACH 

The information-flux approach requires the time- 
evolved form of specific operators Oj in the Heisenberg 
picture, which is given by 

6 ^t) =U{tfd 3 U{t), (2) 

with U{t) the time evolution operator, O = X, Y, Z and 
j = 1,...,N. Using this method, one is able to under- 
stand the dependence of Oj (t) on any Ok and design the 
proper set of interactions in such a way that it becomes 
possible to drive a desired evolution by means of engi- 
neered quantum interference fl2| . In our case, we will 
exploit the time dependence of the couplings in order to 
induce precisely this effect. 

Our task here is the transmission of quantum informa- 
tion from the first spin to the last one in the chain. This 
requires the study of 0jv(i)'s, which can be decomposed 
into the operator basis built out of all possible tensor 
products of { H,. A,. ) ,. Z,\ with 1^ the 2x2 identity op- 
erator applied to spin i. Wc have 

(*o|d*(i)|*o>= E I°°'W(0o|O'il0o), (3) 

0'—X,Y,Z,l 

where |*o) = |0o)i ® ^0)2 jv ^ s the initial state of the 
whole chain, the first spin being in \<po)i and the rest 
of the chain in ^0)2 n (this definition can be easily 
generalized to the case of mixed states of the chain). 
The coefficient I 00 (t) is defined as the information flux 
at time t from 0' x to On [12j |. If our system achieves 
= 1 when O' = and for O = X, Y, Z, we have 
perfect I — > N state transfer. 

We can determine which are the terms of the operator 
basis that are involved in the decomposition of Ojv(t)'s 
by representing such operators in an oriented weighted 
graph, whose construction is straightforward. Each node 
of the graph corresponds to an operator involved in the 
decomposition. If the interest is focused on the evolution 
of On, it should occupy the first node of the graph. Any 
operator resulting from the commutators of Oat with the 
Hamiltonian should be included in the graph and linked 
to it with an oriented edge. These steps are repeated 
as necessary and the construction process ends when no 
new operator is created on commutation. A numerical 
weight, determined by the performed commutations and 
needed in order to calculate the information flux, is at- 
tached to each edge of the graph. Finally, an outgoing 
(incoming) edge corresponds to a + (-) sign. Due to the 
particular form of H in Eq. (Q]), Xn and 3V will gener- 
ate the oriented graph shown in Fig. [T] (where the case 
N = 5 has been considered). Clearly, Yn already appears 
in the graph corresponding to Xn and we do not need to 
construct a dedicated graph for it. 

For the sake of clarity, here we show how to evaluate 
the information flux for a time-independent Hamiltonian 
such as Eq. ([T]) for J x , v {t) — J x . y . The method can be 
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FIG. 1: Oriented graph describing how the operators X5 and 
Y5 evolve in the Heisenberg picture. The operators within 
each circle (a node) give rise to their nearest neighbors on 
commutation with H in Eq. {T}. The oriented edges connect 
such nodes. The corresponding coefficients are also shown 
and an outgoing (incoming) edge with respect to a node im- 
plies a + ( — ) sign. Black, light green and dark red edges are 
associated respectively with B, J x and J y interaction terms. 



When the parameters otj(t) cannot be analytically eval- 
uated due to the difficulties of the evolution, it is still 
possible to approximate them by means of recurrence for- 
mulas. We have 

1=0 

where M is a proper cut-off and 

a 2 } = J y a i~ 1] + Jx(4^ 1] - Bo£j- 1] , (8) 



explicitly extended to any time-dependent model. How- 
ever, in our numerical analysis we divide the total time 
of the evolution in a series of short steps within which 
the Hamiltonian is taken as constant. In this way, in 
the limit of a large number of steps and correspondingly 
short time steps, the flux of information can be accurately 
evaluated by means of only a light computational effort. 
A cut-off in the number of considered temporal steps is 
chosen so as to have no noticeable difference in the re- 
sults if a larger number is taken. For a time-independent 
Hamiltonian, Oj(t) = e int O j e~ lW and, by means of the 
operator expansion formula [15], we have 

4-(*)=Eifer dfc(<5i) (4) 

fc=0 

with Ck(Oj) the nested commutator of order k between 
% and Oj and C (Oj) = Oj. 

Here, for the sake of a clear presentation of our pro- 
posal and without loss of generality, we focus on the case 
of N = 5 (the case of an arbitrary N being easily extrap- 
olated from our example) and study the dynamics of X5 
under the model in Eq. ([1}. The terms in the operator 
expansion are deduced by 

<?i(X 5 ) - [H,X 5 ] = -2iJ y Y 4 Z 5 + 2iBY 5 , 
C 2 (X 5 ) = [H, [H,X 5 ]} = 2iB{2iJ x X A Z b - 2iBX 5 ) 

— 2i J y (2i J y Xc, + 2i J x X 3 Z 4 Z$ — 2i B X 4 Z$), 

(5) 

Clearly, the only operators involved in this iterative se- 
quence are those in the graph of Fig. [TJ Therefore, it is 
possible to write the evolved operators X^(t) as 

X 5 (t) =a 1 {t)X h + a 2 (t)Y 4 Z 5 + a 3 {t)X 3 Z 4 Z 5 + 

a 4 (t)Y 2 Z 3 Z 4 Z 5 + a 5 (t)Xx Z 2 Z 3 Z 4 Z 5 + 

„ , . (6) 
a 6 (t)Y 5 + a 7 (t)X 4 Z 5 + a 8 (t)Y 3 ZiZ 5 + 

a 9 (t)X 2 Z 3 Z 4 Z 5 + a w (t)Y 1 Z 2 Z 3 Z i Z^. 



with otj =0 (1) for j ^ 1 (j = 1). It is immediate to 
note that these recurrence formulas come directly from 
the commutation rules in Eq. ([5]) . They can thus be easily 
derived from the graph in Fig. [TJ The same formulas, 
but with the initial conditions <x- = (1) for j ^ 6 

U = 6), describe the evolution of the operator Y§. The 
approach described so far constitutes the basis for the 
analysis of perfect quantum state transfer via temporally 
controlled information kicking, which is the focus of the 
next Section. 



III. THE IDEAL CASE 

Let us study the action of ti in Eq. (TJJ) when only one 
of the inter-spin coupling terms differs from zero. The 
graph in Fig. [TJ shows that, in these conditions, every 
node is linked to just one nearest neighbor. In fact, the 
full oriented graph for the information flux analysis cor- 
responding to Eq. (TTJ) can be seen as the juxtaposition 
of three mutually disconnected subgraphs, each associ- 
ated with only one of the coupling terms in % and rep- 
resented by a different color in Fig. [TJ It is easy to track 
the operator dynamics associated with these configura- 
tions. Let us consider, for instance, the case B(t) ^ and 
Jx,y{t) = Vt. Fig. [TJ (black color) shows that A5 and Y5 
are mutually linked and, by using the formal quantitative 
analysis of the previous Section, we have 

X B (t) = co S [2p(t)]X 5 + sin[2/?(t)]Y 5j 

(9) 

y 5 (t) = -sm[2p(t)]X 5 + coB[20(i)]y 6 

with j3{t) — J * B(t')dt'. Clearly, an isolated node in a 
graph corresponds to an operator that is constant in time 
[for instance, X 5 (t) = X 5 when J x ^ 0, J y = B = 0, as 
in Fig. [TJ (light green)] . This result is key to our inves- 
tigation. If we let the interaction terms be alternatively 
non-zero for a time window [0,r] such that /3(t) = tt/A, 
the action of % in Eq. (TTJ will correspond to kicks oj 
information between pairs of connected nodes of the 
graphs discussed above. In Fig. [2] we show an exam- 
ple for J x ,y(t) being alternately non-zero for five time 
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intervals such that J t 3 J x (t')dt' — ir/4 for j = 1,3,5 
[ft j _ i J y (t')dt' = tt/4 for j = 2,4], with J x (t) = 

= 0] whenever J y {t) f [J x (t) ^ 0]. The neg- 
ligibility of the mutual overlaps of the coupling functions 
ensures that the information kicks occur according to the 
ordered scheme in Fig. [5] A similar transport of infor- 
mation from the leftmost to the rightmost part of the 
graph is achieved by kicking information using a pattern 
given by fixing J y — (J x — 0) and alternately chang- 
ing the amplitude of J x (J y ) and B in a way so as to 
satisfy the conditions stated above. After a sequence of 
N (for time-dependent J x and J y ) or 2N — 1 kicks (for 
time-dependent B and either J x or J y ), one obtains 

X&(t*) = X1Z2Z3Z4Z5, 

Mr*) = Y 1 Z 2 Z 3 Z 4 Z 5 , ' ' ' 

where r* is the total time of the evolution. By initially 
preparing the j th spin of the chain (except the first) in an 
eigenstate of Zj , there is a complete end-to-end transport 
of information across the chain. This evolution interac- 
tion can thus be exploited for a perfect state transfer ex- 
actly as it happens for the Hamiltonian model of Ref. Q . 
Quite interestingly, any other state-transfer protocol re- 
lying on the same information-flux structure, as those in 
Refs. @, 0], can be mapped into the very same tempo- 
ral pattern of information kicking illustrated here. This 
is the main result of our investigation, proving that an 
alternative to the standard architecture of properly ar- 
ranged interaction-strength patterns for quantum state 
transfer is actually possible. The key is the introduc- 
tion of an alternated series of non-overlapping temporal 
kicks of information. The equivalence of the two scenar- 
ios has here been proven in terms of an isomorphism of 
the respective information fluxes and is, as such, com- 
pletely general. Depending on the details of a specific 
practical realization of the intra-chain couplings, tem- 
poral control of the form highlighted here can be more 
convenient than pre-fabrication of a specific pattern of 
static coupling strengths. This is the case, for instance, 
for quasi-unidimensional optical lattices loaded with neu- 
tral atoms, where inter-site couplings are achievable in 
a time-controlled way via external optical potentials or 
cold atomic collisions. In perspective, this could also 
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FIG. 2: "Operator jumps" induced by H in Eq. l[T]). when J x 
and J y are alternately different from zero. J x (J y ) is non-zero 
in the time intervals ti with an odd (even) value of i, and its 
integral over each of these intervals is equal to 7r/4. 



be a viable option in solid-state structures (such as ar- 
rays of Josephson junctions), where fabrication of exactly 
the pattern required for ideal state transfer would be de- 
manding (if not prohibitive) already at moderate chain 
lengths. Achieving control via proper voltage/magnetic 
pulses inducing time-controlled information kicks would 
be a possibility to exploit instead. 



IV. THE REALISTIC CASE 

In order to go beyond the idealization of the scheme 
assessed in the previous Section and make our analysis 
closer to more realistic situations, the ideal conditions on 
the pattern of time-dependent couplings invoked before 
will be relaxed here. Let us first consider J y (t) = Vt 
with J x (t) and B(t) following the behavior dictated by 

J x (t) = J max [sin(t + n/4)} m , 

B(t) = B max [co S (t + 7r/4)] m , ' ' 

where J m ax = B max are the maximum values of J x (t) 
and B(t) (properly chosen in order to satisfy the integral 
condition) and m is an even number. Pulse shapes of 
this form (for moderate values of m) are routinely gen- 
erated by commercial pulse generators. A large value 
of m obviously implies small overlap between the two 
coupling functions. A period of 2tt will correspond to a 
"horizontal" and a "vertical" kick in the associated oper- 
ator graph. For a numerical analysis of the information- 
flux behavior under the action of this time-dependent 
Hamiltonian, we divided the total time interval [0, 2Nw] 
in short steps. Within each of them, both J x (t) and B(t) 
are taken as constant and numerically equal to the av- 
erage values they assume in the respective time window. 
In this way, the recurrence formulas presented in Sec. |TT] 
allow us to estimate the evolution of Xn or Yn for the 
duration of the step. Clearly, the shorter is the time step, 
the more accurate is the estimate. The method is then 
iterated for each step and the total evolution of Xn or 
Yn is finally reconstructed. As the number of operators 
involved in each decomposition is equal to 2N, a light 
computational power is required. 

In order to evaluate the performance of our scheme, 
both the coefficient ajv(t), attached to X1Z2 ■ ■ ■ Zn in 
the decomposition of A?jv(t), and ci2N(t), attached to 
Y1Z2 ■ ■ ■ Zn in the decomposition of $Ar(£)j should be 
studied. However, one can equally well characterize the 
whole state transfer process simply by considering cvjv (t) 
only. The maximum value of ajv(i) for two values of m is 
studied against N in Figs. [3] (a) and (b). It is matter of 
performing an explicit calculation to see that the maxi- 
mum of a2N{t) follows similar behaviors. We then quan- 
tify the mean transmission fidelity F(N, t) by averaging 
the state fidelity (ipm\PN{t)\ipm) between a pure input 
state IV'in) to transfer and the state Pat(<) of spin iV at 
time t. By assuming a uniform distribution of \tpi n )'s [3[, 
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FIG. 3: Maximum value of Qjv(t) against the number of spins 
N in the chain. The value of m in Eq. is respectively 
equal to 4 [panel (a)] and 6 [panel (b)]. 



we find 
F{N 1 2Nn)= l - 



I + n N (2.V7r) ( | + iajv(2iV7r) 



(12) 

For instance, if we take m = 6, we obtain for any length 
of the chain N < 25 a value of F(N, 2Nn) > 0.984. 
As ajv(2./V7r) — > 1, the average fidelity converges to 1, 
resulting in perfect state transfer. 

The ability to temporally control both J x and B can 
however be experimentally demanding. In particular, 
while the arrangement of external potentials with a de- 
sired time-dependence can be relatively straightforward, 
this might not be the case for J x (t)'s, which are in general 
determined by coupling mechanisms internal to the chain 
itself and hardly controllable with the degree of accu- 
racy required by our scheme. We thus make our requests 
lighter by taking J x as constant and allowing a time de- 
pendence only for the external magnetic field. Obviously, 
it is formally the same to take a constant B with time- 
varying J x (t). However, pragmatically our choice seems 
to us to be the least demanding. We thus adapt the 
strategy assessed before to the case of J x — J cons t- Our 
quantitative study shows that a slight modification to 
the sinusoidal trend described before is sufficient to re- 
tain the efficiency of the protocol. We take well-spaced 
square-shaped pulses for B(t) as shown in Fig. 01 The pa- 
rameter (5 quantifies the width of magnetic field pulses. 
It can be roughly seen as half of the overall period di- 
vided by the time in which B is close to its maximum 
value. Clearly, a large value of 6 corresponds to a more 
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FIG. 4: Temporal behavior of B (red line) and J x (blue line), 
for a value of 5 = 5. The quantities plotted in the vertical 
axis are in units of J. 
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FIG. 5: Maximum value of a/v(t) against the parameter 5, 
for a chain with N = 5 spins. 
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FIG. 6: Maximum of ajv(t) against the parameter 8, for N = 
2k + 1 with k = 2 — ► 7 (going from the top to the bottom 



ideal case, as the overlap between the two functions is 
smaller, hence affecting the maximum value that ajv(f) 
can achieve and, in turn, the transfer fidelity. The depen- 
dence of aN^Nir) on 5 is shown in Fig.[5]for N — 5. The 
fidelity F(5, IOit) is already larger than 0.99 for 6 = 8. 
Of course, the effect of the overlap between the two func- 
tions depends on N: the longer the chain, the lower the 
transmission fidelity. The dependence of the fidelity on 
the length of the chain is plotted in Fig. [51 We note 
that, for iV up to 15, a value of S = 16 is large enough 
to obtain the maximum of ajq ~ 0.96 (corresponding to 
F ~ 0.974). 

To give a complete overview of the protocol, one can 
finally analyze, for a fixed 5, the behavior of the max- 
imum of (x n against N. The case of 6 = 20 has been 
considered, and the results are shown in Fig. [71 In these 
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FIG. 7: Maximum of arf(t) against the length of the chain TV 
for 5 = 20. 
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conditions, the corresponding transmission fidelity for a 
chain with N — 25 spins is ~ 0.947. 

We have also studied the case in which the non-zero in- 
teraction terms are J x and J y , with B = 0. The transmis- 
sion fidelity is larger than the one obtained with non-zero 
Jx (Jy) and B, as the number of kicks required is smaller 
and the effect of the overlap is thus reduced. However, 
we presented the results for non-zero J x (J y ) and B as 
we consider the choice of a time-dependent magnetic field 
the least demanding from an experimental point of view. 



V. GENERATION OF MULTI-PARTITE 
ENTANGLEMENT 

In this Section we discuss another interesting aspect 
of the proposed model, namely the possibility to gen- 
erate genuine multi-partite entanglement. Bi-partite as 
well as multi-partite entanglement creation in spin chains 
has been recently studied from different points of view. 
For instance, one can generate entanglement by means 
of a proper pattern of coupling parameters [l4[ , bydriv- 
ing the evolution through a resonant interaction [16| . by 
performing a global quench [17j , or by leaving two empty 
sites in a uniformly filled chain 18]. By assuming the 
ability to initialize the state of the chain before the ap- 
plication of the pulse sequence in the protocol proposed 
in this article (we necessarily have more demanding re- 
quirements for this particular task), Greenberger-Horne- 
Zeilinger (GHZ) states can be generated Such a 

resource is well-known to be useful for multi-agent pro- 
tocols for distributed QIP like quantum secret sharing, 
remote implementation of unknown operations and quan- 
tum average estimation [20|. Let us consider again the 
Hamiltonian % in Eq. ([1} under the ideal conditions men- 
tioned above. It is interesting to notice that, after the 
proper kicking of information, the evolution of all the 
spin operators associated with the elements of the chain 
have exactly the same form as those obtained by using 
the model proposed in Refs. [!, Q. 

We now want to evaluate the dynamics of a particular 
pure state of the chain under the action of this Hamil- 
tonian. We will need to abandon the assumption that 
at most a single spin excitation populates the system. It 
is straightforward to check that, whenever spin i is in 
an eigenstate of Z^, the corresponding symmetric spin 
N — i + 1 will end up in the same state. Whenever two or 
more spins are initially in eigenstates of X or Y , the cor- 
responding symmetric spins will end up in an entangled 
state. Let us consider, for definiteness, an even number of 
spins. We call {di, 0,2, a n } the indices corresponding 
to spins initially prepared in eigenstates |d) a . = |0, 1) 

of Z and {61, 62, b m } those labeling spins initially pre- 



pared in eigenstates \d)bi = |+7~ )&. of X (the analysis 
can be straightforwardly generalized to the case of odd 
number of spins and eigenstates of Y). The initial state 
of the chain is thus 



|*o}= (g) \d,d) ai , bi 

{ai ,bi} 



(i = 1, ...,n). 



(13) 



It is matter of a straightforward calculation to see that, in 
this case, the final state corresponds to a mirror-inversion 
operation on the state 



!*(*)) = 



1 

71 



{<*<} 



{b,} 



bi +*(-!)' ^)Z bi \d) t 

{b,} 



(14) 

with Z = 0, 1 depending on the actual form of |^o)- 
Clearly, \^(t)) is the tensor product of a separable state 
for spins {a\, a 2 , a n } and an m-particle GHZ-like state 
for spins {61, b 2 , b m }. 



VI. REMARKS 

By means of the information flux approach, we have 
presented a new strategy for obtaining perfect state 
transfer in a finite, open chain of spins. We have shown 
that the evolution of spin operators under the action of 
a site-dependent Hamiltonian that is already known to 
allow perfect state transfer can be mimicked by the use 
of a homogeneous time-dependent Hamiltonian. The en- 
gineering of the coupling strengths, necessary in the first 
scenario, is thus converted to a temporal arrangement of 
the interaction terms. We have discussed a particularly 
interesting case where such a time dependence can be 
reduced to the requirement of just a non-constant mag- 
netic field affecting the spins of the chain. Interestingly, 
this very same model can also be used in order to create 
a GHZ-like entangled state shared by a sub-set of spins 
of our choice, leaving all the other spins in a separable 
state. We believe that our proposal to achieve perfect 
state transfer and quantum correlation sharing in many- 
body registers of interacting spins widens the possibili- 
ties for efficient and non-demanding short/medium-haul 
quantum communication. 
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